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Abstract 

The time evolution of a charged spin 1 /2 particle interacting with a magnetic field 
is analyzed in the framework of the complex spectral theory, based on Prigogine's 
principles that provide a rigorous description of irreversibility. A detailed survey 
of the irreversible relaxation process of the spin 1/2 particle is carried out on the 
basis of the subdynamics theory. We obtain the results for the Markovian and the 
non-Markovian evolution of the charged particle. 
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1 Introduction 

Our experience suggests that the symmetry in the time is disrupted and the 
future and the past play different roles. The world surrounding us has obvi- 
ous irreversible nature. However, Poincare showed that [1] the description of 
the irreversibility is impossible on the basis of the classical laws of dynamics, 
since the latter are reversible in time. The equations of quantum mechanics 
are reversible too. The description of the physical world on the basis of fun- 
damental classical and quantum theories as defined by the laws of the nature 
is determined to be time reversible. Difference between the classical descrip- 
tion of the nature and those processes in the nature which we observe creates 
the conflict situation. Therefore the problem of the description of irreversible 
world on the bases of the reversible equations of classical and quantum physics 
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is raised. I will not discuss the solutions of this problem in classical physics, I 
will examine the quantum mechanics. First of all it is necessary to note that 
the question about irreversibility of quantum processes can be examined in 
connection with the problem of interpretation of quantum mechanics. Briefly, 
I will examine two most important approaches - Copenhagen and Everett's 
interpretations. 

Copenhagen interpretation of quantum mechanics asserts that we cannot speak 
about the quantum properties of the system before these properties are mea- 
sured [2], PJ. In other words, quantum theory describes not quantum world 
but the fact that we can speak about quantum world after measurement. Let 
the energy of the quantum system has two values Ex, E 2 and <px, ¥2 corre- 
sponding wave functions. If the system is described by the wave function ip 
it is possible to present it in the form of the superposition of the functions 
tpx and ifi2'- if > = a ifi + 0'2 l f2- Before measurement the wave function ip oc- 
cupies simultaneously two levels and energy of the system does not have a 
specific value. Only after the measurement we will obtain the values E\ or E2 
with probabilities | ai | 2 , | a 2 | 2 respectively. The passage from the "poten- 
tial possibilities" which are described by the wave function ip to the "actual 
realities", which can be measured, is called reduction or collapse of the wave 
function. The physical process leading to collapse is called quantum decoher- 
ence. The process of the measurement leads to the abrupt, irreversible change 
of the state. Any measurement is not a reversed process. Thus, irreversibil- 
ity appears when the measurement is carried out. In this case the behavior 
of the quantum system is described by Neumann's postulate of reduction jl] 
which, however, does not answer the question which dynamic processes lead to 
irreversibility. Thus, in accordance with Copenhagen interpretation of quan- 
tum mechanics the observed processes proceeding in the world are caused by 
our measurements including, for example, the quantum transitions. However, 
as Prigogine notes, since the quantum transitions are the basic mechanism of 
chemical reactions it is difficult to agree with the latter assertion. Is it possible 
to consider the chemistry as the result of our observation? If yes, then who 
observed the chemical reactions which led to the appearance of life [5]? 
Another interpretation of quantum mechanics is proposed by H. Everett (rela- 
tive state interpretation) [H] . It was developed in the works [7] , [H] (many- world 
interpretation); [9], [TO] . [TT] (extension of Everett's concept). Many-world in- 
terpretation allows the existence of the infinite numbers of classical realities: 
each term of quantum superposition corresponds to one classical world. Ac- 
cording to the extended Everett's concept the measurement entangles the 
measured system with the environment, however, linear superposition is not 
destroyed. Thus, after measurement the linearity of quantum mechanics is pre- 
served. In this case the reduction of the wave function does not occur. How- 
ever, at the specific moment of time the observer sees only one result of the 
measurement - the one classical world. This occurs because the consciousness 
of the observer divides the state of quantum world into the classical realities 
which he receives independently. The consciousness of the observer subjec- 
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tively makes the selection of one alternative. The consciousness continuously 
creates classical reality [10J. Hence follows: the classical world is the illusion, 
reduction of the wave function is the illusion too, since they appear subjec- 
tively in the consciousness of the observer. In the approach, the concepts of 
the "arrow of time" and irreversibility are rather connected with its subjective 
feeling of the observer then with the objective property of material itself. It 
is assumed, the further study of this interpretation will be possible after the 
construction of the model of " quantum consciousness" , which is discussed in 
the work [TT] . 

Obviously, the approaches examined above contain a subjective, anthropo- 
morphic element. The interpretations contain a basic distinction between the 
quantum system and the observer. The latter fact is undesirable for the theory 
since it takes us away from the solution of the basic problem: the determina- 
tion of the objective laws of the nature (if we believe in the objectivity of 
the nature). If Copenhagen interpretation is the direct consequence of the in- 
sufficient realization of the physical sense of the basic postulates of quantum 
mechanics, then Everett's interpretation is a natural consequence of them. 
This situation impels one to the alternative formulation of quantum dynamics 
which does not appeal to the observer excluding the subjective, anthropomor- 
phous element. In the approach the irreversibility must be represented as the 
property of material itself and is not defined by the active role of the observer. 
In the paper I examine the alternative formulation of quantum dynamics in 
the framework of the Brussels-Austin group works that has been headed by 
I. Prigogine for many years. In the works, a possible variant of description 
of nonequilibrium processes at microscopic level in the frame of the Liouville 
space extension of quantum mechanics is investigated. The mechanism of the 
asymmetry of processes in the time which made it possible to accomplish a 
passage from the reversible evolution to irreversible one was developed. Thus, 
new irreversible dynamics with the disrupted time symmetry was formulated. 
The symmetry in the time disrupted as a result of asymmetric nature of the 
physically permissible states. 

At the present moment, it is necessary to continue further development of 
the Brussels - Austin group approach in the framework of the realistic models 
of interaction. Then, the irreversible evolution of a charged spin 1/2 particle 
interacting with a magnetic field is investigated on the basis of subdynamics 
theory. In the work the survey of the basic principles of the Brussels-Austin 
group is carried out. In section 2 the eigenvalue problem is discussed. Per- 
turbative solution of the Schroedinger equation in the framework of complex 
spectral representation is given. In section 3 the Liouville formulation of quan- 
tum mechanics is represented. The task of the complex spectral representation 
of the Liouvillian is solved in section 4. In section 5 the theory of subdynam- 
ics is discussed. The time evolution of the density matrix is determined. The 
interacting model of positive charged spin 1/2 particle with magnetic field is 
examined in section 6. The complex eigenvalue problem for the spin system is 
solved in section 7. The expression for the density matrix element describing 
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the irreversible evolution of the charged particle spin system and numerical 
calculations are obtained in section 8. 



2 Perturbative solution of the Schroedinger equation, complex spec- 
tral representation 



As the first step, in the framework of Schroedinger equation, I examine as 
the asymmetry in time can arise. I examine the eigenvalue problem for the 
Hamiltonian H = Hq + XV 

H | 1p a >= E a I V> a >, (1) 

where H Q - free Hamiltonian, XV - interaction part, A - coupling constant. 
In the conventional case, Hamiltonian if is a Hermitian operator, E a is a 
perturbed energy of the state - a real number. It is known that the usual 
procedure of equation ([!]) solution on the basis of perturbation method can 
lead to the appearance of the denominators l/(E a — E a >), where E a , E a i 
are the energies corresponding to the unperturbed situation. Obviously, the 
divergences can arise at E a = E a i. The basic question now is - what we can 
do to avoid the divergences, when E a = E a >. I examine the situation using the 
simple Friedrichs model [12] (the model is presented closely to the text of the 
works [13] - [IB]). Despite the fact, that the solution of the problem for the 
Friedrichs model is known it serves as a good example for the demonstration 
of the essence of situation. The model describes interaction of two level atom 
and electromagnetic field. In the Friedrichs model | 1 > corresponds to the 
atom in its bare exited level [H], | k > corresponds to the bare field mode 
with the atom in its ground state. The state | 1 > is coupled to the state | k > 

H = H + XV 

= Ex | 1 >< 1 | + J2 E k \kxk \ +XY,Vk(\ k >< 1 | + | 1 >< k |), ( 2 ) 

k k 

where 

| 1 >< 1 | +^2 | k >< k \— 1, < a | a >= 8 oaf. (3) 
k 

Here a (a') = 1 or k. In the ordinary approach the eigenvalue problem for the 
Hamiltonian H is formulated as follows 

H | V>i >= Ei | V»i > ■ (4) 
For the eigenstate | ipi > (f° r small A) perturbation method gives the expres- 



sion 



^>«|i>-Ei^Vl*>- ( 5 ) 
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If Ei > 0, divergences appear at E^ = E\. I. Prigogine and co-workers noted 
that in the general case the satisfactory solution of this problem is impossible 
on the basis of the conventional formulation of quantum dynamics. However, 
the eigenvalue problem can be solved if the time ordering of the eigenstates will 
be introduced. This procedure can be realized through the introduction into 
the denominators imaginary terms: — ie for the relaxation processes, which 
are oriented into the future and +ie for the excitation processes, which are 
oriented into the past. In this case the eigenvalue problem Q is reduced to 
the complex eigenvalue problem 

H\<p 1 >=Z 1 \ ipi >, <(p x \H =< (pi | Zi, (6) 

where we must distinguish right - eigenstate | <p\ > and left - eigenstate 
< (pi | [oTll3|ll4] . Zi is a complex 

Zi=Ei- i 1 . (7) 

Here E% is a renormalized energy and 7 is a real positive value. This pro- 
cedure makes it possible to avoid the divergences and leads to the following 
expressions for the eigenstates | <fi >, < (pi | [13J 

\<Pi >«|1>~E jr, X / k , + I k >, (8) 

k \ h k - hi- Z)_ h 

<yil«<i|-E fF f k z)+ < k \- 0) 

k \ h k - hi - Z)_ h 

In the expressions ^ and ^ the designation l/(Ek — E\ — z)\^ has been 
referred to as "delayed analytic continuation" [T3]TH] and can be determined 
through the integration with a test function f(Ek). For example, 



{E k - Ei - z)t iry z^-i-y\J Ek-Ei - zJ z£C +' 

where we first have to evaluate the integration on the upper half-plane C + 
and then the limit of z — >• —27 must be taken. Thus, the approach makes it 
possible to solve the problems which could not be solved within the frame- 
work of classical and quantum mechanics. Now we can realize the program 
of Heisenberg - to solve task at the eigenvalues which could not be solved 
within the framework of traditional methods. The spectral representation of 
the Hamiltonian takes the form 

H = ^2Z a \ (p a >< (p a 1 . (11) 

a 

For the eigenstates we have relations 

J2 I Pa >< ¥a |= 1, < \ W >= <W- (12) 
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Since, H is Hermitian the corresponding eigenstates | <px >, < 0x | are outside 
Hilbert space, have no Hilbert norm 

< <fx | ipx > = < 0x | 0i >= (13) 

and furthermore have a zero average energy 

< (px | H | (ft >=< 0x | H | 0x >= 0. (14) 

The eigenstates | tpx >, < 0i | ar e called "Gamow vectors". Extensive litera- 
ture is dedicated to the study of " Gamow vectors" , for example [T7] - [20] . 
Thus, the Hermiticity of H leads to the fact that "usual" norms of eigenstates 
| <fx >, < 01 | disappear. However, the eigenstates | <px >, < 0i \ have a bro- 
ken time symmetry. We can associate | (fx > with the unstable state, which 
vanishes for t — > +oo, | 0x > corresponds to the state, which vanishes for 
t — > — oo 

| ipx(t) >= exp(-tExt - yt) | ^i(O) >, (15) 
| 0x(t) >= exp(-iExt + -ft) | 0x(O) > . (16) 
However, the are few unpleasant aspects. 

Using the eigenvectors | <px >, \ 0i > it is possible to construct the density 
operators [H] 

Pa =| 01 >< (fl |, Pb =| (Pi >< 01 | ■ (17) 

They are invariants of motion and they do not describe the process of relax- 



ation. In addition, "unusual" norms and zero average energy (13), (14) exist. 
These are the reasons why a more general space of density matrix-Liouville 
space will be examined. Nevertheless, it is important to emphasize: the com- 
plex perturbative solution of the eigenvalue problem in the framework of the 
Schroedinger equation makes it possible to determine 7 which has important 
physical sense. The value 27 determines the probability of transition between 
two quantum states [H] , [21] . This result will be important for us in the next 
sections. 



3 Liouville formalism 



Now I will examine the equation which determines the density matrix p - 
Liouville-von Neumann equation 

i% = L P . (18) 

Liouville-von Neumann operator (Liouvillian) has the form: 

L = H x 1 - 1 x H, (19) 
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here symbol " x" denotes the operation (A x B)p=ApB. Operator L is written 
down in the sum of free part L Q , that depends on the free Hamiltonian H and 
the interaction part Lj, that depends on AV^ = Hj - Hamiltonian of interaction: 
L = Lq + Lj. Let state \a > be a eigenvector of the free Hamiltonian, then 
we have the equation H \a >=u a \a > with the energy ui a . The dyad of the 
states | a >< (3\ is the eigenstate of operator L 



L \a >< (3\ = {u a - ujp)\a >< (3\ (20) 

or 

L \a;P)) = w a>/3 \a;(3)), (21) 

where the designations \a >< (3\ = \a; (3)) and w a ^ = u) a — up were used. 
In the Liouville space for the dyadic operators we have the relations: 
the inner product defined by 

({A\B))=Tr(A^B), (22) 

the matrix elements are given by 

<a\A\(3>= «a;/3L4», (23) 

the biorthogonality and bicompleteness relations have the form: 

((a;/3'|a;^))=W W , X>;/3)>«ar;/3| = l. (24) 

For operator L we have the equation (the text is written close to the materials 
of the works [5]fT4"] ) 

L | f a ,p)) = W a ,(3 | fa,(3)), (25) 

here | f a ,p)) =| 4> a >< ^8 | ; ^a,/3 = E a — Ep - the eigenvector and eigen- 
value of Liouville-von Neumann operator L, where | ip a > (| tpp >), E a (Ep) 
- the eigenvector and eigenvalue of the total Hamiltonian ([T]). The spectral 
representation of operator L has the form 

L = Y,\ fa,(3))WaA(f<*,f3 I • (26) 

Since w a , a = E a — E a = 0, the diagonal elements of density matrix do not 
depend on the time. Therefore, strictly speaking, the conventional formalism 
does not make it possible to describe the evolution of the system to the equi- 
librium. Eq. (25) can be rewritten by using the correlation index v. 

L\U))=w v \U)), (27) 



where w u = w a $ ( or w u = w a ^ for the case (21)) and | f v )) =| f a ,i3))- 
The details of the theory of correlations can be found, for example, in the 
works [5,22.23]. In equation (27) i/ = 0ifa = /3 - that is the set of diagonal 
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operators | a >< a | - "vacuum of correlations" and v ^ is the remaining off- 
diagonal case. Then, the concept of degree of correlation must be introduced. 
The degree of correlation of dyads | u)), d, has been defined as the minimum 
number of interactions Lj by which a dyadic state can reach the vacuum of 
correlation. For example, the degree of vacuum of correlation d = 0. The con- 
cept of correlation will be examined lately in connection with the model which 
is investigated in the article. 

Certainly, the problem of the description of irreversibility can be solved on 
the basis of well known physical approximations [24J. The solution is based 
on the distinction between the open system and its environment. The envi- 
ronment is assumed to be in the thermodynamic equilibrium. However, as it 
was noted [25] "This distinction introduces an anthropomorphic element. In- 
deed, the environment, introduced by one observer, may be different to that 
seen by another observer. Moreover, this distinction is phenomenological as 
one avoids answering the most fundamental question of nonequilibrium sta- 
tistical mechanics, namely, what is the mechanism that governs the approach 
to thermal equilibrium of the environment, without violating the basic laws 
of physics." The approach of Brussels-Austin group is based on the assertion: 
the irreversibility is the objective property of quantum world. This objective 
reality we have to describe without including the anthropomorphic principle, 
such as lack of our knowledge, coarse-graining approximations, environmental 
approach [26]. The solution of the problems, which was examined above, is 
proposed on the basis of the Liouville space extension of quantum mechan- 
ics [2], [IB] or otherwise, on the basis of the complex spectral representation 
of Liouville-von Neumann operator L, and subdynamics theory. In this case, 
Liouvillian has complex eigenvalues, that break time symmetry without intro- 
duction of an anthropomorphic element or dissipative terms to the Hamilto- 
nian [27] or to the Liouvillian |28j . 



4 Complex spectral representation of Liouvillian 

We have a new formulation of the eigenvalue problem: 

L | = Z ] | | L = | Z;, (28) 

where Z v - are the complex values, v is a correlation index and j is a degener- 
acy index, since one type of correlation index can correspond to the different 
states (the complex eigenvalue problem for the Liouvillian is examined in the 
works [14,29,30]). It was shown [T4l23|l31] that the eigenstates of L can be 
written in the terms of kinetic operators C u and D v . Operator C u creates cor- 
relations other than the v correlations, D u is destruction operator. The use of 
the kinetic operators allows to write down the expressions for the eigenstates 
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of Liouville operator in the following form [T4] 



= W) 1/a ^ol«i». = ((^|^(iV;) 1/2 , (29) 

where 



% = P v + = + D u , (30) 

here iVJ - is a normalization constant. The determination of the states 

and operators P u , C u , D u can be found in the works |T4 .15. 22J . In the 
general case, for example, the operators P u are determined by the following 
expressions [22] 

p v = £ K»«^l, (W)> = <wv (si) 

Let me examine some results obtained in the original works. The work [T4] 
shows that § U D satisfy so-called nonlinear Lippmann-Schwinger equation. 
For the $^ we have 

$^ = F+VF- — 1 [L/^ - $ u r Li$ u r }P u , (32) 

where the time ordering —ie^ u was introduced. This is the important moment 
of the theory because the time ordering breaks time symmetry and, thus, leads 
to the description of the irreversibility. The determination of the sign of the 
infinitesimals is connected with degree of correlation d. In accordance with 
our experience we will consider the directions to the higher degrees of corre- 
lation oriented in the future, for example, the decay of elementary particles, 
quantum transitions with emission or absorption of energy, and the directions 
to the lowest degrees of correlation are oriented in the past. This leads to the 
relations: 



= +£ if df, >d u (t> 0); e^ u = -e if < d u (t < 0). (33) 
For the we have the equation 

1 



= p" + P V [^ V D L! - & v D Li$ v D ] Yl 



W v — — IE 



(34) 



I'/i 



D v as follows 



Eqs. (32), (34) determine the kinetic operators of creation C v and destruction 

[Lj&c-^Lj^P", 



c u = J2 p 



-i 



D v = P V W D U - ® V D L& V D ] £ P^ — 



(35) 



(36) 
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It is evident from the last expressions that D v ^ {C v )\ but D v = (C u )*, 
where "f" corresponds to the Hermitian conjugation and the action "*" - the 
"star" conjugation which is Hermitian conjugation plus the change e 



Substituting the expression (29) in Eq. (28) and multiplying P u from left on 
both sides, we obtain [H] 



Oc\<))=Z»\u»)), 



where 



P U L(P U + C u ) = L P U + P v Li{P v + C U )P V 



(37) 



(3J 



Also, it is possible to show the validity of the following useful relationship 



(39) 



The latter will be necessary for us in the future. In Eq. (37) Q v c is the collision 



operator connected with the kinetic operator C v . This is non- Hermitian dis- 
sipative operator which plays the main role in the nonequilibrium dynamics. 
As it was shown in the work [22] operator 9q can be reduced to the collision 
operator in Pauli master equation for the weakly coupled systems. Comparing 



Eqs. (28), (37) we can see that \u u ,)) is eigenstate of collision operators Q v c with 



the same eigenvalues as L. 

Thus, the time ordering introduces into the theory the "arrow of the time", 
and leads to the connection of quantum mechanics with kinetic, time irre- 
versible dynamics. As an example, I examine the Friedrichs model. In the 
works |14j . [T6] the following results were obtained: 



L | *S» = Z\ | ¥?», 
Z? = -2i7, 1*0)) = 



[v 



0), 

yi;yi)>; 



L | *")) = z u | * y )), (u^O) 

ylj/c y Tp 7k,l TP 

L — L\ — H/k, Zj — r^k — 



7k±' 



where, for example, 



Eu — Eu , 



^ k ' k ')) =| (p k ;(pv)). 



(40) 



(41) 



(42) 



In (40)-(42) c.c. is the operation of the complex conjugation, Z\ is determined 
by relationship Q, | (fx > is determined by (J8]) and | tpk >, | <pk > can be 
found in the works [2], [TB] and correspond to the field mode. 
Obviously, the complex eigenvalues Zj of the Liouville operator cannot be 
reduced to the simple difference of two eigenvalues of the Hamiltonian as in the 

are not the multiplications of 



case (25). Furthermore, eigenstates 



the wave functions, they cannot be obtained from the Schroedinger equation. 
In this sense we have the alternative formulation of quantum theory. 
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5 The theory of subdynamics, time evolution of the density matrix 



Equation (18) does not describe the irreversible evolution. The precise de- 
scription of irreversibility requires the determination of the new approach. 
This approach was developed on the basis of the theory of subdynamics. In 
the theory of Brussels-Austin group, subdynamics is called the construction 
of a complete set of spectral projectors IT [T3~lll4|ll6j (see also [32|33|34j ) 



n" = EI*i»«*il- 

3 

The projectors IF satisfy the following relations 

WL = LW, (commutativity); *^W = 1, (completeness); 
WW' = W5 UU >, (orthogonality); W = (IT)* , (star-Hermiticity) 
and in our approximation 

W — > P u if interaction is absent (A=0). 
Operator W can be represented in the following form 

W = (P v + C U )A U (P U + D U ), 
where A v is the star-Hermitian operator A v = (A u )* 

A u = pu^p V + d u c v )- 1 p u . 



(43) 



(44) 



(45) 



(46) 



(47) 



Taking (44) it is possible to write down the density matrix p in the form 

pW = En>(t) = E^), 



(4£ 



where p u = II^p. Using (jl Sp , (]39|), (|44|) and (|46|) it is possible to write down 

(49) 



p v (t) = Wexp(-iLt)p(0) = 
[P v + C v ) exp{-t$ u c t)A"{P u + D u )p(0). 



I examine the component P u p v (t). It was called the "privileged" component 
of p u (t). For the P u p u {t) it is easy to obtain 

P u p v {t) = P v e^{-i^ c t)A u {P v + D u )p(0), (50) 

where in our case p(0) = P°p(0) - vacuum of correlations and the relations 
P V C V = 0, P U P V = P v were used. The relation (50) leads to the equation 

.dP v p v {t) 



dt 



r c p» P v {t). 



(51) 
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Thus, it is shown operators p v satisfy separate equations of motion. The equa- 
tion (51) is the kinetic equation for each U u subspace. This is the reason why 
the projectors IF can be associated with the introduction of the concept of 
subdynamics. 

Using (48) I determine the component P°p(t) 

P°p(t) = P°J2 P V {t) = P°P°(t) + P° E P V (t). (52) 



The expression P°p (t) is determined by (50) for v = 0. The second term has 
the form 



p ° E P V ( t ) = £ P °( pU + ° U ) exp(-t^ u c t)A u (P u + D u ) x 
p°p(0) = Y,P°C"exp(-ir c t)A' / D u P p(0). 



(53) 



The obtained expressions lead 



P°p(t) = P° exp{-i#° c t)A o P°p{0) + E P ° C " exp(-i^t)A v D v P°p(0). ( 54 ) 

Taking the time derivative of P°p{t) - component we have 
.dP°p(t) 



Of 



^?7(i)+ 



J2 P°C u r c exp{-ir c t)A u D u P°p{0). 



(55) 



The equation (55 ) consists of two parts: the first one corresponds to the Marko- 
vian approximation, the second part includes the memory effects and deter- 
mines the non-Markovian processes [25] . 

It is possible to conclude, the introduction of the subdynamics makes it possi- 
ble to determine the Eqs. (51) and (55) which are the fundamental equations 
for the description of the irreversible processes. 



6 Charged particle spin behavior in the magnetic field 



I examine the behavior of the positive charged particle in the magnetic field [35] 
As the charged particle it is possible to examine proton, positron and any pos- 
itive charged particle with non-zero magnetic moment and with spin 1/2. 
Spin S of the particle is connected with magnetic moment ft by the expression 

P = 1 P S, (56) 
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where the values j p , S are determined by 

7 P = ^, S=\a. (57) 

Here e, m p are the charge and the mass of particle, g p is ^-factor, a - Pauli 
matrices. I use the system of units, where h = 1 and the speed of light c = 1. 
Let the isolated charged particle with magnetic moment /x be placed into the 
external magnetostatic field which is directed along z axis - H z . The Hamil- 
tonian operator then becomes 



p 



In the expression (58) a z is a z-component of the vector a. 

Particle with spin 1/2 in the external magnetic field H z has two energy levels 

corresponding to two different values of magnetic quantum number m 

E m = —m-fpHz, where m = ±1/2. (59) 

The distance between the levels is determined by the energy AE 

AE = u = lp H z . (60) 

The transitions of the particle from one energy level to another E_y 2 ^ £1/2 
can occur as a result of the influence on the system of the external magnetic 
field Hi perpendicular to the field H z . The exciting magnetic field Hi must be 
variable and its frequency uj must coincide with frequency ojq. Such transitions 
are accompanied by emission or absorption of energy AE. This is the well 
known phenomenon of the magnetic resonance. 

Strictly speaking, the energy levels are not determined accurately. They have 
an uncertainty of the values. It leads to the spectrum of the radiated energies, 
which distribution is simulated by the introduction of distribution function 
f(oo) with the width A. The distribution function describes Lorentzian or 
Gaussian distribution of energy. Lorentzian distribution is described by the 
function [35] 

/M = ; L. r! 7V2 > 6 = l A ( 61 ) 

7T A + {UJ — UJq) z I 



with norm 



f{uj)duj = 1. (62) 



Taking into account the distribution f(u) it is possible to determine the rate 
of the transition between the energy levels [35], [36] 

Wi/2^-1/2 = ^-1/2^1/2 = W= l -^lHlf{u). (63) 
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In the simplest case, the components of the field Hi change according to the 
law [37] 

H x = Hi cos(W), H y = — Hi sm(ojt), (64) 

where u is the frequency of the rotating field in the xy plane. Then, the 
Hamiltonian operator takes the form 



H 



magnetic 

e 

^ P 4m„ 



-9 P 



Arrir 



a-H 



C X H X + UyHy + <7 Z H Z ) . 



(65) 



Whence it follows 



H 



magnetic 



-9 P 



H z H T — iH„ 



H x + iH v 



H 7 



Since, 

H x — iH y = Hi exp(iujt), H x + iH y = Hi exp(-iujt) 
I obtain the expression 



(66) 



(67) 



H 



magnetic 



Hsh, + H, 



'9 P J— a zH z - g p 







Hi exp(iu)t) 



Arrir 



Am p \Hiexp{-iut) 



(6h 



SH X 



The first term Hsh z m the expression (68) determines interaction with the 



field H z , the second term Hsh xv determines interaction with the field Hi. 
I examine the field Hi so that Hi -C H z (see A. Abragam [35J). This condition 
has the greatest practical interest. In this case for the eigenfunctions m , which 
are determined by the equation 



Hmagnetic^m ^mtymi 



m = ±1/2, 



(69) 



I take the approximation 



h/2 = | | exp(-iEt/ 2 t), 



-1/2 



exp(-iE-y 2 t). 



(70) 



Then, the function, which describes the behavior of spin, can be determined 
by the superposition of the functions <px/2, 0-1/2 



O1/201/2 + O-1/20-1/2 



Ol/2 



1 

I exp(-iEi /2 t) + a_i /2 I I eacp(-i£7_i/ 2 t). 



(71) 
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For the Hermitian conjugated (h.c.) function ip h ' c ' I have 



</2 (10) exp(iE 1/2 t) + a h _:t/2 (01) exp(i£_ 1/2 t) 



(72) 



The expression (68) corresponds to the monochromatic case. For the passage 
to the general case it is necessary to make the following replacements 



Hi exp(ioot) 



A(u)Hi exp(iut)dtu, 



(73) 



Hi exp(— iut) — > J B(Lu)Hiexp(—icot)dto, 

— oo 

where the value of the coefficients A(u), B(uS) will be determined later. 



7 Complex eigenvalue problem for the spin system. 



I examine the transition E-1/2 — > ^1/2 • This transition is accompanied by 



the emission of the energy AE (60). The transition with the emission of the 



energy is the irreversible process. The description of the irreversible evolution 



of the spin system will be carried out on the basis of the expression (54). I 
will determine the evolution of matrix element < — 1/2 | p(t) \ — 1/2 > which 
corresponds to the probability of finding the particle in the state with spin 
— 1/2 at the time moment t. 

The evolution of the system in the kinetic - Markovian approximation has 
the exponential nature ^ exp(— wt). So, in the work [21] (in this work the 
irreversible evolution of the unstable ir~ - meson was defined) w = 2 / ~f Pw - the 
rate of the ir~ - meson decay. I assume that in our case w = 2 / -f_i/ 2 = W, where 



W is the rate (63) under the resonance condition lu = lu (lu is determined 
by (60)). Therefore, the first step will consist of the determination of the 
value 7-i/2- This task will be executed on the basis of perturbative solution 
of complex eigenvalue problem in the framework of the method of the second 
quantization. 



Let me examine the replacement (73). According to the method of the second 



quantization I define the values A(oj), B(oS) as operators 

A(oo) - g(u)at(u), B(u) - g(u)a(u). (74) 

Here aJ(u), a(u) are the operators of creation and destruction of photon with 
the energy u. I assume, these operators satisfy the commutation relation 

[a(u),a\u')) = 6(u-u'). (75) 
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In the determination (74), g(u) is the function satisfying the condition 



g 2 {uj)duj = 1. 



(76) 



As will be shown below (see (101)) the square of the function g(tu) is the 
distribution function 



and 4> h ' c ' in (72) we must make a±u 2 



The coefficients a±i/ 2 in the expression (|71| find 

h.c. 
±1/2 

are 



the value of the destruction operators of particles with m = ±1/2. For a hA 



h.c 



a 



ih.c. 



±1/2' 



here a 



the creation operators. It is assumed, the operators a±i/2, 
anti-commutation relation 



a 



t 

±1/2 



±1/2 

satisfy the 



[ a m, <4n'] + 



6 r , 



m (m') = ±1/2. 



(77) 



In this case the Hamiltonian of the particle-magnetic field system can be 
determined in the form 



H 



magnetic 



4m T 



$<T X H X 



( 



fjp 



4m r 



OO \ 

/ g(u})aJ(u})Hi exp(iut)duj ' 



/ g(u)a(u)Hiexp(—iut)du 

\— oo 



(78) 



where the previous designations of operators are preserved. Now, the values 



(f>, (jr and Hamiltonian (78) can be examined as the operators in the Dirac 
representation. 

I examine the eigenvalue problem for the Hamiltonian H = Ho ± Hj, where 
a free part H = ± Hsh z , and an interaction part Hj = XV = Hsu xy - 
is a free electromagnetic Hamiltonian. Let state | uo >= a*(u) | > (| > is 
the vacuum state: a{uS) \ >= 0) be a eigenvector of the free Hamiltonian H u 
then, \ uj >= uj \ uj >, where u> is the photon energy. The Hamiltonian H w 
can be determined as follows 



ujcr (uj)a(u)du. 



(79) 



As will be shown below (appendix C), for our model, here and in the sub- 
sequent expressions integration over u have to be carried out in the interval 



to £ [0, oo). I determine the states | ±1/2 >= 
are the eigenvectors of the Hamiltonian Hsh z 



*±l/2 



>. The states | ±1/2 > 



H 



SH Z 



-9 P 



4m r 



E 

m=l/2,-l/2 



JO) 
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and Hsh z | ±1/2 >= E±i/ 2 | ±1/2 >. In this case the expression H \ 
±1/2, cu >= (u + E±i/2) | ±1/2, to > is valid, where | ±1/2, to >= a} ±1 , 2 a\ijo) \ 

>. Since, the transition (E_y 2 — > E 1 / 2 ) occurs as the result of interaction 
with the field if i that is why the term H$ h x is defined as the interaction part 
XV. 

1 will solve the problem assuming, that the eigenvalue Z_i/ 2 of the Hamil- 
tonian H is complex (the general formalism of the complex spectral repre- 
sentation can be found in the works [T3~|14lll6|l38f39| ). In accordance with the 
approach [13], in our case, I will distinguish equation for the right eigenstate 
| (f-i/2 > arid for the left eigenstate < ip-1/2 | of Hamiltonian H 

H I y?-i/2 >= ZL1/2 I ^-1/2 >, < I H =< ip-1/2 I Z-i/2- (81) 

I expand the values | ^-1/2 >? < ^-1/2 I, ^-1/2 in the perturbation series 



V-1/2 



>= E A n l A >, < $-1/2 1= E A n < tf% /3 1, 



(n) 



n=0 



n=0 



^-1/2 = E ^"^-l/V 



52) 



n=0 



where 



-1/2 >, < <p% 2 \=< -1/2 



7(0) 

Z -l/2 



E-1/2 



A = e — the charge of particle. 



34) 



In accordance with the definitions (81 )-(84) for the first equation of (81 ) I can 
write down 



< -1/2 I (H J2e n I <p% 3 > +eV E e" | tp% 3 > 



(n) 



n=0 



n=0 



(JXJ LXJ 

<-l/2|Ee"^ 2 Ee"'|^ 2 >, 



35) 



where | — 1/2 > is our initial state. 

Hence, it is easy to obtain the expression for the coefficient Z_^, 2 



r(n) 
'-1/2 



< -1/2 I V I > - E ^}/2 < -1/2 I > 



n-l 



1=1 



„(n-Q 
1/2 



36) 



term in ( 78 ) leads to 



Let us examine the first term in the expression (86). The use of the interaction 

9 P Hi 



< -1/2 1 v 1 ^ >= 



4m 



- f g(u) exp (-i(E 1/2 - E_ 1/2 + u)t)x 

n J 



(87) 



< 1/2, uj I V 9 ^!/^ > duj. 



/2 
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The substitution of the expression (87) into (86) and summing up Yl e " on 

n=l 

the both sides of the obtained expression lead to the result 



OO TT °9 OO 

1^% = -^ / sHE^expH(£i /2 -£-i /2 + ^) 

n=l 4m P _4 n=l 

oo n— 1 

< 1/2, w | ^ >da;-Ee fl E Z<? 1/a < -1/2 | > 

n=l i=l 



X 



whence I obtain 



Z-i/2 — E-i/2 



q H °r 



< 1/2, w | ^-1/2 > rfcj - (Z_i /2 - S_i /2 )(< -1/2 | ^_i/ 2 > -I] 



J9) 



Let us determine the value < 1/2, u | <f-x/2 >. For this purpose I write down 
the obvious equality 



< 1/2, uj I (H £ e" I <^ /2 > +eV £ e» | <p% > 



1/2 ^ I V-1/2 

n=0 n=0 



<l/2,^|Ee^i n 1 ) /2 E^I^2> ) 



n=0 



n'=0 



(90) 



where | 1/2, u > is our final state. 



The expression (90) leads to 



< 1/2, u I 



•1 



Ex/2 +UJ — E_i/ 2 



f < 1/2,0; 1 v|^T/2 ) >-E^!/2< 1 /2 



I (n-0 ^ 
I V-l/2 > 



(91) 



It is assumed, the state with the spin —1/2 disappears in the future, converting 
to the state with the spin 1/2 and photon with the energy u. In this case the 



time ordering of the expression (91) has to be introduced. For this purpose I 
make the replacement 



1 



1 



Ex/2 + U — E-x/2 E X /2 + 10 - E-i/ 



2 - i£ 



(92) 



Here the symbol \i in infinitesimal e^ u corresponds to the state | 1/2, uj >, the 
symbol v corresponds to the state | —1/2 >. The state | 1/2, uj > has the high 
degree of correlation and I assume e^ u = e > 0. 
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For the first term of (91) I have the result 



< 1/2, u; | V | <p%$ >= -|^(u;)exp(i(£ 1/2 -£_ 1/2 + u/)t)x 



< -1/2 | yfc? > . 



(93) 



The substitution of the expressions (93), (92) into (91 ) and summing up Y< e " 

n=l 

lead to 



n J-f 

< 1/2, u | <p-i/2 >= e^—^g(u) exp(i(Ei/ 2 - E^x/2 + v)t) x 

1 

< -1/2 | <p_ 1/2 > 



(94) 



E1/2 + to — Z_y 2 — it 



Taking into account the result (94), the expression (89) can be represented in 
the form 



Using the formal expression 
1 



1 



E1/2 + u - Z_i/2 - is 



-duj. 



(95) 



Ex/2 + CO — Z-i/2 — IE Ex/2 + U — E-i/2 

I present the eigenvalue Z_i/2 as follows 

Z-X/2 = #-1/2 — H-l/2i 

where E_x/2 is a renormalized energy 



iir5(Ex/2 + uj - E_x/2) (96) 



(97) 



1 



^/-^-(^) 2 p7,v)^ +w _ b _ i/2 

P - the principal part of the integral 



-du. 



and 



7-1/2 = ^l p Hx9 M- 



(99) 



Comparing the obtained expression (99) for j-x/2 with the expression (63) for 
the rate W we see 



W-X/2-*l/2/u 

where the equality 



27-i/2/ff3(wo)=/(«o) = i^lpHxfM^ 



(100) 
(101) 
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is used. Evidently, in this case, the obtained expression 27_!/ 2 coincides with 
the expression for W_i/2-»i/2 under the resonance condition u = P_i/ 2 — Pi/2- 



8 Charged particle spin evolution 



Let us examine the matrix element ((-1/2; -1/2 | P°p(t))) (see the expression 
(54)). The projection operator P° (31 ) can be redetermined in the form (A.l ). 
For the matrix element the following designations are valid 



((-1/2; -1/2 | P<°>p(f)» = ((-1/2; -1/2 | p(t)» 
< -1/2 | p{t) | -1/2 > . 



(102) 



For the operator A v (47) the equality 



A v = P U (P U + D u C u Y l P u = P V + Y / (-l) n (D u C u ) r 



(103) 



n=l 



is correct. Then, in the first approximation, it is possible to assume that 

A v « P u . (104) 



Using the expressions (54), (102) and (104) I obtain 

((-1/2; -1/2 | p(t))) = ((-1/2; -1/2 | exp(-<t)P°p(0))) + 
£((-1/2; -1/2 I CeM-^ct)P v D v P°p(0))). 



(105) 



First of all, I examine the first term in the expression (105). This term cor- 
responds to the Markovian approximation. In our case the initial condition is 
given by 



P(0) 



4/2 >< -1/2 |, i.e., ((-1/2; -1/2 | p(0)» = 1, 



and for our matrix element I have 

((-l/2;-l/2|exp(-<t)P°p(0))) = 

((-1/2; -1/2 | exp(-z^) | -1/2; -1/2)). 



From the expression (38) it follows 



here it was assumed that 



•&% = P LjC°P , 



P L!P = 0. 



(106) 
(107) 

(108) 
(109) 
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According to [32] operator C v (the expression (35)), in the first approximation 
can be represented in the form 



y P m — 



is 



110) 



Hence, for the operator C° I obtain 



LjP 



(111) 



Here the symbol v — corresponds to the vacuum of correlations, p in in- 
finitesimal e^o corresponds to the high degree of correlation. Infinitesimal 
determines the time ordered transition from the vacuum of correlations to 



the high degree of correlation. In our model the initial condition (106) cuts 



from the possible vacuum of correlation states the initial correlation state 
| —1/2; —1/2)) which corresponds to initial pure state | —1/2 >. High degree 
of correlations /i are determined by the correlation final state | 1/2 u; —1/2)) 
which corresponds to the final state | 1/2, a; >. In this case = e > 0. The 



details of the evaluation of the expression for the matrix element (107) are 
examined in the appendix A: 



((-1/2; -1/2 | exp(-^)P°p(0))) = exp(-2 7 _ 1/2 £). 



112) 



The non-Markovian effect is determined by the second term in the expres- 



sion (105). For it I have the result 



£((-1/2; -1/2 | C v eM-Wct)P V D v P p(V))) 



eg p H x 
4m n 



,exp(-z(o; -uj )t), 

f(u)—f ^^ duJ + c - c - 

(u> — lo — is) 1 



(113) 



The details of calculations of the expression (113) are represented in the ap- 
pendix B. 



Taking into account the determination (102) and the results (112), (113) the 



time evolution of the charged particle spin system can be represented by the 
expression 



< -1/2 | p(t) | -1/2 >= 



exp(-27_i /2 t) + 



egpHi 
4m„ 



exp(-i(w - uo)t) 
j(u)— , Nn doj + c.c. 



( u — ujq — xe ) 



114) 
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The calculation of the integral in the expression (114) (the features of the 
calculation are examined in the appendix C) leads to the final result. 



< -1/2 | p(t) | -1/2 >= 
exp(-2 7 _ 1/2 t) + 



A(t) sm{uj t) + B(t) cos(uj t) ) . 



115) 



The coefficients A(t), B(t) in (115) are determined by the expressions 

a 4 (a 4 (l + b) - a 2 (l + 66)£ 2 + 6£ 4 ) 



A(t)= /deexp(-0- 



oo 

B(t)= |d£exp(-£) 



D(0 

-a 7 (2 + 46)£ + 4a 5 6£ 3 



The function D(£) has the form 



D(0 = a 8 (l + b) 2 + a 6 (2 + 6(2 + 46))^ 2 + a 4 (l + 6(-2 + 66))£ 4 + 



2 f 8 



116) 



117) 



118) 



a 2 6(-2 + Ab)C + &T 
where a = Uot and 6 = (u /5) 2 . 

Thus, the privileged component in the determination (52) corresponds to the 



kinetic, Markovian process (the first term in (115)), the non-Markovian effect 



coming from the nonprivileged component (the second term in (115)) 



Figure 1 presents the calculations of the separate contributions: the Markovian 
and the non-Markovian terms. Calculations are executed for proton: proton 
mass - 938.27 MeV, g-factor - 5.58, nuclear magneton - 3.15 x 10~ 14 MeV T _1 
[36] , [40J. It is seen, the relaxation process is determined by the Markovian 
term. The kinetic evolution prevails in the entire time interval. The non- 
Markovian term is determined by the damped oscillations and is not essential. 
The other sets of the values (H z , Hi, 5) substantially decrease the contribution 
of the latter. Note, that the calculation of the non-Markovian effect was car- 
ried out earlier, for example, in the work [26] (a quantum harmonic oscillator 
linearly coupled to a bosonic massless scalar field), where the contribution of 
the latter is more noticeable (see also [2S], [II])- A time scale of transition of 
the evolution from the non-Markovian regime to the Markovian regime was 
determined. Thus, it was shown, the non-Markovian term corresponds to the 
well-known quantum Zeno effect [42]. 



9 Conclusion 



In conclusion it is necessary to note that the approach I have examined is not 
the one more additional interpretation of quantum mechanics. It is not in- 
tended to interpret the quantum mechanics, but it is proposed to reformulate 
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it so that the irreversibility naturally would enter into quantum dynamics. 
The approach denies the conventional opinion that the irreversibility appears 
only at the macroscopic level, while the microscopic level can be described 
by the laws, reversed in the time. Thus, the irreversible nature of quantum 
physics is asserted. In this case, the basic problem, designated by Boltzmann 
and Planck - to formulate the second law of thermodynamics at the micro- 
scopic level, is revealed. In the approach the solution of this problem is pro- 
posed in the works [2], [30], [S] - [46J. The law of the increase of entropy is 
accepted as fundamental there, that determines the "arrow of time", the dif- 
ference between the past and the future. The future corresponds to the larger 
value of the entropy. Thus, the existence of the "arrow of time", the second 
law of thermodynamics at the microscopic level are determined as the basic 
postulates. Unconditionally, the interaction with the environment leads to the 
irreversible effects, but in the Brussels-Austin group approach the environ- 
ment is not limited by the special requirements as, for example, to be in the 
thermodynamic equilibrium. The preceding calculation of the relaxation time 
(Fig.l, Markovian term) for a spin 1/2 particle interacting with a magnetic 
field demonstrates the utility of the Brussels-Austin formalism in a practical 
problem that may be tested in experiment. 

Is it possible to speak about the laws of quantum world independently of the 
observer? Bohr came to the conclusion about the impossibility of description 
of the independent quantum processes. The approach of the Brussels-Austin 
group is the attempt to learn to describe the processes in the quantum world 
as the property of the material itself without active role of the observer. 
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Appendix A. Markovian approximation 

In our approximation the projection operator P° (31) {y — - vacuum of 
correlations) can be redetermined in the form 



P°=|-l/2;-l/2»«-l/2;-l/2| + 
1 



1 (A 



where the volume f2 is determined so that ^ = 1 and 6(0) is the delta- 
function. For the operators P v [y ^ - high degree of correlations ) I determine 

oo 

P v = I I -1/2; l/2w))((-l/2; l/2w | du+ 

— oo 

oo 

J I l/2u; -1/2))((1/2cj; -1/2 | dcu. 



(A.2) 
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Then, the following relationships are valid: 



PT...P U = P°, P V P V ...P V = P u , P K) P V = P V P U = 0. (A.3) 



The expansion of the expression (107) in the series leads to 



((-1/2; -1/2 | exp(-i^t) | -1/2; -1/2)) = 
1 + (-it)((-l/2; -1/2 | $° c | -1/2; -1/2)) + 

^£((-1/2; -1/2 | (^) 2 | -1/2; -1/2))+ ( A - 4 ) 
^£((-1/2; -1/2 | (^) 3 | -1/2; -1/2)) + ... . 



Taking into account ( JlOS] ) for the term ((-1/2; -1/2 | $° c | -1/2; -1/2)) I 
obtain 

((-1/2; -1/2 | ^ | -LA -1/2)) = 
((-1/2; -1/2 | P%C°P | -1/2; -1/2)) = 

((-1/2; -1/2 | U £ P" ~\ Lj | -1/2; -1/2)) = (A . 5 ) 
((-1/2; -1/2 | Li £ P** ~\. L t 1-1/2; -1/2)), 

where e^ = e > 0. In our case 

pa,p -1 = 

l/2u- -l/2»«l/2c; -1/2 | ^du+ (A Q) 

Wl/2w,-l/2 — V ' 7 

-l/2;l/2w))«-l/2;l/2w | — Aj. 

W-l/2,l/2u> - ^ 

Here, w 1 / 2lu -i/2 = — W-i/2,i/2w = ^1/2 + a; - P_i/ 2 - Using the determinations 
L 7 = ffj x 1 - 1 x H T , (Ax B)C = ACB; (A.7) 

((a; (3 I U I a'; /?')) =< a | ifj | a' > Spp - <W <P'\H I \p>, (A.8) 

where 5^ (5 aQ /) is delta-function if the indices /?', /? (a, a') correspond to the 
continuous spectrum and it is Kronecker's symbol if the indices correspond to 



-00 
00 
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the discrete spectrum, for (A. 5) we will have 
((-1/2; -1/2 | 0° c | -1/2; -1/2}} = 

oo 

- J < -1/2 | Hi I 1/2, w >< 1/2, uj \ H T | -1/2 > 

— oo 

\du). 

Wl/2u> -1/2 + 



Wl/2w,-l/2 - 



(A.9) 



Keeping in mind the result (97) it is possible to determine the equality 

1 1 



El/2 +LO — Z-l/2 



IE 



1 



(E 1/2 + u - z)+ i7 _ i/2 
0(e), 



(A. 10) 



E1/2 + uj — .EL 1/2 — 
where 0(e) is the term of higher order on e. 

In the equality (A. 10) the concept of "delayed analytic continuation" is used 
(section 2). From the relationship (A. 10) another equality follows 

1 1 



E 1/2 + u-Z±f /2 + ie (E 1/2 + u-z) 

+ 0(e). 



+*7-l/2 



(A.ll) 



Ei/2 + uj - E-x/2 + is 



In this case the operation of integration of the expressions, which contain the 



values of the form 



(E 1/2 +u>-z) +i , 



is determined by the rule: we first have to 



-1/2 



evaluate the integration on the lower half-plane and then the limit z +27-1/2 
must be taken. Taking into account (A.10), (A.ll ), using the Hamiltonian (78) 
and limiting by order e 2 it is easy to obtain 

((-1/2; -1/2 I 0° I -1/2; -1/2)} = 
feg p H 



fegpHiY 
\ 4m„ / 



00 


9 2 


— 00 




00 




h 2 




—00 





(E 1/2 + u-z)± 



-du)+ 



(E 1/2 + 



du. 



(A.12) 



UJ — z 



+*7-l/2 



Comparing the obtained expression (A.12) with the expression (95) and taking 
into account the designations (97), (A. 10), (A.ll) I can write down 



((-1/2; -1/2 I V% I -1/2; -1/2)) 

Z-l/2 — Z-X/2 = —^1-1/2- 



(A. 13) 
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Making the analogous calculations for the each term of the expression (A.4) 
in the limit Q — > oo I will obtain the result (112). 

Appendix B. Non-Markovian term 

In this appendix I examine the details of the obtaining expression (113). The 
determination of the operator C v (110) leads to 

(C) 



W v — — IE 



(B.l) 



where the action " *" corresponds to the "star" conjugation (section 4). Taking 



ment (101), in the lower-order approximation for the operator exp(— i^t) 
(see (J38J)) 

exp(-ztf&*) « exp{-iw l/ t)P u (B.2) 



it is easy to obtain 

£«-l/2; -1/2 | ^exp(-z^)P^P°p(0))) = 

((-1/2; -1/2 | LjP u exp(-M u c t)P u Lj | -1/2; -1/2)) 



E 

WexpH(£„-£_ 1/2 )t) 
p ^ 

exp(-z(£_i /2 - E p )t) 



[W — IE) 



< -1/2 \ H I \p><p\H I \ -1/2 > 



(E p - £_ 1/2 - z £ ) 2 
< p | #/ | -1/2 >< -1/2 \Hi\p> 

(E-l/2 



- E p - ie) 



(B.3) 



y (exp(-i(E 1/2 + u - E-!/ 2 )t)x 

— oo 

< -1/2 I Hi I 1/2, u >< 1/2, u \ Hi | -1/2 > 



(E 1/2 + uj- J5_i/a - ie) 



da; + c.c. 



4m„ 



/(a;) — r^^duj + c.c. ) . 



[U — Uq — IE) 



Here, index p corresponds to the state | 1/2, u >, E p = E 1 / 2 + u, u = 
.EL 1/2 — Ei/ 2 . In (B.3) the sum corresponds to summation (integration) over 
all discrete (continuous) indices. 

Appendix C. Calculation of the integral 



The calculation of the integral in the expression (114) will be carried out in 
the kinematic region, where the function f(u) satisfies the condition 

f(u)du = 1. (C.l) 
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In this case, the integration can be realized in the interval u G [0, oo). 
Since, t > 0, we have to deform the contour of integration in the lower half- 
plane [23]. In this region the application of the Cauchy's Theorem is possible 
(see Figure 2) 



J exp(—izt)F(z)dz = J exp(—iut)F(u)dw+ 
r R o 

o 

J exp(-izt) F \z)dz + / exp(— izt)F(z)dz = 0, 



1R 



-iR 



where 



F{z) 



[z-u; o r 



If f(z) is the Lorentzian distribution (61) then 



lim / exp(— izt)F(z)dz — ► 0. 



It leads to 



(C.2) 



(C.3) 



(C.4) 



exp ( — iuit ) F (u ) dui 



exp(—izt)F(z)dz. 



(C.5) 



The replacement z = —iy, where y — £/t makes it possible to reduce our 
integral to the form 



exp(-i(u - u )t) 
lu — luq — isy 



. . . exp(-i(-iy)t) 

-« / dyfy-W) , .... — —r^ exp(iuj t) 



-iy - u ) 



(C.6) 



.exp(^ot) , 

- z — / 



exp(-0 



l + 6(l + if) (l + if) 



where a = wot and 6 = (cg>o/5) 2 . It allows us to evalute the obtained integral 
with the use of the " Mathematica" -program. This leads to the expression 



j{Ld)— — h c.c. I doj 



[UJ — LOq — IE) 



(C.7) 



A(t) sin(a;ot) + B(t) cos(uj t) 



The coefficients A(t), B(t) in (C.7) are determined by the expressions (116), (117) 
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Fig. 2. 

Figure Legends 

Fig.l: the Markovian and the non-Markovian terms; the proton, H z = 10 4 Oe, 
E x = 10 2 Oe, 5 = 0.033 m _1 . 

Fig.2: a contour rotation. 
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